In-plane optical features of the underdoped La 2 Cu04 based compounds: Theoretical 

multiband analysis 
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The three-component a&-plane optical conductivity of the high-T c cuprates is derived using the 
gauge invariant response theory, and compared to the data previously obtained from the optical 
reflectivity measurements in the La2Cu04 based families. The valence electrons are described by 
the Emery three-band model with the antiferromagnetic correlations represented by an effective 
single-particle potential. In the < S < 0.3 doping range, it is shown that the total spectral 
weight of the three-band model is shared between the intra- and interband channels nearly in equal 
proportions. At optimum doping, the low-frequency conductivity has a (non-Drude) nearly single- 
component form, which transforms with decreasing doping into a two-component structure. The 
mid-infrared spectral weight is found to be extremely sensitive to the symmetry of the effective 
single-particle potential, as well as to the doping level. The gauge invariant form of the static and 
elastic Raman vertices is determined, allowing explicit verification of the effective mass theorem and 
the related conductivity sum rules. 
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I. INTRODUCTION 

Following the marginal Fermi liquid theory pj, the 
nested Fermi liquid theory or the ordinary memory- 
function approach Q, the low- frequency excitations of 
the electronic system involved in the optical conduc- 
tivity and electronic Raman scattering spectra of the 
overdoped high-T c cuprates are described in terms of 
the frequency-dependent effective mass m(u>) and the 
frequency-dependent relaxation time t(u>), with a mi- 
nor role played by selection rules or by normal-state 
coherence factors. Attempts to explain the measure- 
ments in the underdoped cuprates using similar single- 
component models have been unsuccessful. Anomalous 
low-frequency structures related to the superconducting 
(SC) pseudogap, the antiferromagnetic (AF) pseudogap, 
the two-magnon excitations or the density waves are ob- 
served in the latter crystals, usually in one or two (of 
four) "optical" channels (optical conductivity and A-\ „ , 
Big and B 2g Raman symmetries) . 0, 0, IE M IE S E3 
Any theory of the low-frequency excitations in the under- 
doped cuprates should explain anomalies in both optical 
conductivity and Raman scattering spectra, taking the 
normal-state coherence effects properly into account. 

The previous theoretical work on the optical conduc 
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tivity of the underdoped cuprates 
is carried out for the one-band models with a particular 
care devoted to the treatment of strong local correlations 
responsible for the AF structure in the low-frequency con- 
ductivity, but no attention was paid to the (interband) 
optical processes in the visible part of the spectra. The 
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principal physical problem is that when such one-band 
models are applied to the Raman spectra, then the reso- 
nant nature of the Raman scattering processes is ignored, 
because the intermediate states in the Raman scatte ring 
processes are treated in the static approximation. [17lll8| 
The first experimental evidence for the crucial role of the 
interband processes at fkv 1.75 — 2.75 eV for the low- 
frequency physics of the high-T c cuprates was given by 
Cooper and coworkers Q , considering the correlation be- 
tween the optical excitations across the charge-transfer 
gap in YBa2Cu307_j; and the resonant behaviour of the 
magnetic peak in the B\ g Raman response. Similar (but 
less pronounced) dependence on the frequency of the in- 
cident /scattered photons is also found in the Drude part 
of the Raman spectra. The latter is also accompanied by 
several anomalies regarding the relative spectral weights 
of three Raman polarizations. [1 U ^} 

Actually, the recent interest in the Drude part of the 
optical conductivity and Raman scattering spectra of the 
underdoped compounds is stimulated by the clear evi- 
dence for important role of the normal-state coherence 
effects. The optical conductivity and B 2g Raman spec- 
tra scan the electronic states in the nodal region of the 
Fermi surface (so-called "cold spots" on the Fermi sur- 
face: k x — k y ), while the A\ g and B\ g Raman scattering 
processes probe the electronic states in the vicinity of the 
van Hove points ("hot spots"). A striking difference be- 
tween these two groups of spectra occurs for the nearly 
half filled conduction band in all cases where there is 
a perturbation peaked at (usually commensurate) wave 
vector Q close to the nesting vector of the Fermi surface 
2kp. The AF fluctuations at small dopings and the 
charge-density-wave (CDW) fluctuations at the 1/8 dop- 
ing 8] are typical examples. There is also a surprising 
similarity between the mid-infrared (MIR) peak in the 
optical conductivity, /ia>p 0ak , and the magnetic peak in 
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the Big Raman spectra. Both appear at nearly the same 
energy, exhibit similar doping dependence and possibly 
have the same physical origin. 0,0, El El 

It is interesting also to note that in the materials where 
the MIR peak appears at a relatively small energy, com- 
parable with the typical damping energy E=S/t0,||, 
there is an essential difference between the predictions 
of the usual one-band optical models and measurements. 
It is explicitly shown that the development of the mea- 
sured MIR structure with temperature can be fitted well 
with the orbital Kondo model of Emery and Kivelson 
[T^ | (representing the regime where hjujp 0ak « S) rather 
than with the mentioned one-band models (underdamped 
regime: liWp Cak is proportional to characteristic correla- 
tion energy scales, such as the Hubbard interaction U 
[Hill or the AF exchange energy J 

The main goal of the present analysis is the following: 

(i) to determine for the Emery three-band model |2jj all 
electron-photon coupling functions (various current and 
Raman vertices, to be defined below) important for the 
understanding of most of the aforementioned anomalies, 

(ii) to show that, due to the gauge invariance, the present 
multi-component optical model allows a natural transfor- 
mation of the spectra from the underdamped regime into 
the overdamped regime, |2l| and (iii) to point out that 
the correlation effects peaked at Q « 2kp (the AF cor- 
relations in the cuprates, or the umklapp scattering pro- 
cesses in the Bechgaard salts 0, H3, 0) affect the 
optical conductivity of two-dimensional (2D) and quasi- 
one-dimensional (Q1D) systems in a similar way, leading 
in both cases to the dramatic decrease of the effective 
number of conduction electrons. 

However, to ensure a simple analytic form of the cur- 
rent and Raman vertices, the present analysis will be 
limited by the use of several approximations, (i) The 
effects of the AF correlations on the response functions 
will be included in an adiabatic way. (ii) The strong lo- 
cal correlations associated with the Hubbard term on the 
copper sites will be treated in the mean-field approxima- 
tion (MFA) , with the same (renormalized) copper-oxygen 
transfer integral describing the low- and high-frequency 
excitations, (iii) As usual for the in-plane features of 
the high-T c cuprates, the 2D representation of the elec- 
tronic system is used, neglecting several important ques- 
tions dependent on the electron propagation in the di- 
rection perpendicular to the conduction planes, such as 
the anomalous c-axis conductivity or the dimensionality 
crossover. 

In this paper, the response of the valence electrons in 
the high-T c cuprates to the external vector potential is 
analyzed by means of the Emery three-band model [20| . 
The AF correlations in the nearly half filled conduction 
band are modeled by an effective single-particle inter- 
action A(k) peaked at the commensurate wave vector 
Q = (-7T / 'a, tv I 'a). In the calculation we adopt a gauge 
invariant approach which has been proved to enable a 
systematic study of the multi-component optical conduc- 
tivity spectra |25| and which is a simple generalization 



of the memory- function approach [2l|. It will be shown 
that an anisotropic-s effective interaction results in two 
distinct features in the low-frequency optical conductiv- 
ity, with the relation between these two features and the 
high-frequency excitations across the pd charge-transfer 
gap which quantitatively agrees with the experimental 
observations (including the position of thresholds and the 
ratio between the respective spectral weights). 

In Section 2 we describe the simplest version of the 
three-band model which focuses on the optical excita- 
tions across the pd charge-transfer gap and on the low- 
frequency excitations across the MIR correlation gap 
2A(k) (in the rest of the text the MIR conductivity re- 
lated to A(k) will be referred to as the MIR gap struc- 
ture), and review other work. We determine the gauge 
invariant form of the static and elastic Raman vertices 
(Appendices A and B), which makes possible the con- 
sistent formulation of both the three-component optical 
conductivity and the related Raman correlation func- 
tions. Section 3.1 discusses the conductivity sum rules 
of the present three-component optical model, while in 
Section 3.2 the dependence of the low-frequency spectra 
on the symmetry and magnitude of the effective inter- 
action A(k) is studied. Finally, the mutually competing 
effects of the intracell hybridizations and the strong local 
correlations on the interband spectra of the lightly-doped 
conducting planes are briefly discussed in Section 3.3. 



II. THEORETICAL MODEL 

The electrodynamics of valence electrons in the high-T c 
cuprates is described in terms of the extended Hubbard 
model (i.e. the Emery three-band model [2(j in which all 
Coulomb interactions are neglected except for the Hub- 
bard interaction on the copper sites Ud)- As mentioned 
in Introduction, various one-band versions of this model 
have been widely used to describe qualitatively the low- 
frequency optical conductivity. Let us first briefly review 
these models. Then follow the detailed explanation of 
the present optical model and the comparison with the 
experimental data. 

The theoretical study of the underdoped regime reveals 
the MIR gap structure in the optical conductivity, in ac- 
cordance with experimental observation. In the one-band 
Hubbard model jllL Il2l ] the corresponding threshold en- 
ergy is proportional to the Hubbard interaction U (pre- 
cisely, ^Wp Cak ~ U/2, with an additional, charge-transfer 

gap at fou™cak ~ ^Oj wnne m the t-J model this 
energy is related to the AF exchange energy J. Most 
of the relevant results are calculated using the numeri- 
cal diagonalization procedure at high temperatures (usu- 
ally T > 500 K). Similarly, the anomalous (non-Drude, 
single-component) optical conductivity measured in the 
overdoped compounds can be explained in terms of the 
marginal Fermi liquid theory [lj or the nested Fermi liq- 
uid theory |2j, which are characterized by strong quasi- 
particle damping effects. It is essential to notice that fits 
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of the related generalized Drude formula to the measured 
spectra require usually large damping energies (typically 
n/r(0.5 eV)« 0.5 eV 0,01). 

In order to treat both the local correlations related to 
Ud and the interplay between various small energy scales 
more accurately (at low temperatures, with presumably 
small damping energies), as well as to describe the rel- 
evant intracell processes explicitly, different three-band 
versions of the Emery model are examined and contrasted 
to the measured low- and high-fre quency optical spectra 
and the related Raman spectra [l^lMl^lMl^l30ll3l| . 
Kotliar, Lee and Read [26| have calculated the ground 
state energy of the electronic system self-consistently us- 
ing the slave-boson procedure. They have shown that 
the characteristic of the low-frequency optical processes 
is the renormalization of the copper-oxygen bond energy 
t p d (i.e. the copper-oxygen transfer integral, in the usual 
language) . This renormalized bond energy is found to be 
strongly dependent on doping, in particular in the lightly 
doped region, where t p d oc y/S (hereafter, 6 is the hole 
doping). In the case where J = 0, the MIR structure is 
found to be related to the renormalized splitting between 
the oxygen and copper site energies . ( i.e. the renor- 
malized pd dimerization potential) |26J, |2jj. However, in 
the complete slave-boson model (27J this scale can be as- 
sociated with the AF exchange energy J as well. In both 
cases, the excitations across the charge-transfer gap are 
connected with the bare pd dimerization potential. 

The electrodynamics of the underdoped cuprates may 
be close to that of the Bechgaard salts (TMTSF) 2 X 
(X = PF 6 , C10 4 ). Namely, for the nearly half filled 2D 
or Q1D conduction band, the SC and AF (or spin-density 
wave) regions in the phase diagram are expected to be 
controlled by small changes in the already small effective 
number of conduction electrons caused by the interplay 
between the corresponding correlation energy and vari- 
ous other small energy scales. [2l| The measurements 
on the high quality crystals |22L l23| have shown that 
the electrodynamic features of the Bechgaard salts are 
strongly affected by the umklapp scattering processes, in 
particular by the competition between the characteris- 
tic umklapp energy and the transfer integral in the di- 
rection perpendicular to the highly conducting direction 
[23L |24| . In the present paper it will be proposed that 
the two-component low-frequency optical conductivity of 
the underdoped cuprates is connected with a similar cor- 
relation energy related to the AF pseudogap processes, 
and that the transformation of this two-component form 
into a single-component form, in the overdoped regime, 
is related with the interplay between this scale and the 
energy difference between the Fermi energy and the van 
Hove energy. By analogy with the Bechgaard salts, one 
should extend the present analysis by considering the 
competition between the characteristic AF energy and 
the oxygen-oxygen transfer integral t pp (this is expected 
to be relevant for the in-plane conductivity in the lightly- 
doped regime), or the interplane transfer integral t± (im- 
portant for the c-axis conductivity and the dimensional- 



ity crossover problem); however, these are beyond the 
scope of the present work. 

Finally, it is important to note that the present optical 
conductivity analysis, together with the related Raman 
analysis [l8| , can be easily connected with the Hall coeffi- 
cient measurements 1321 , due particularly to the analytic 
form of the intra- and interband current and Raman ver- 
tex functions. It should also be noticed that the model 
parameters extracted below from the optical data can be 
confirmed by angle-resolved photoemission spectroscopy 
(ARPES) experiments |M H HI IM El H| (the con- 
duction band is wide and nearly half filled) , as well as by 
electric field gradient (EFG) measurements [3^, |4jj (the 
copper-oxygen hybridization is relatively strong, again, 
supporting the picture of wide bands |41j)- 

A. Bare three-band Hamiltonian 

The present response theory is based on an effective 
single-particle description of the valence electrons. The 
analysis starts with the 2D three-band Hamiltonian of 
the form [HIH 

H = y^Hg, 

L 

H ° = E^( k )4^k.. (i) 

In the hole picture, the bonding band (the band index 
L = D) is nearly half filled, while the antibonding and 
nonbonding bands (L — P and N) are empty. 

Two versions of the model will be considered. The 
model A includes only two parameters, the splitting be- 
tween the oxygen 2p CT and copper Zd x i_ y i site ener- 
gies A p d and the average first-neighbor bond energy t p d- 
Since two different limits of the three-band model, the 
Ud = limit and the MFA of the Ud — ► 00 limit (here- 
after, the large Ud case), are represented by this model, 
the parameters A p d and t p d can be associated either with 
the bare parameters or with the parameters renormalized 
by the large Ud (in the latter case, the present model is 
identical to the original slave-boson model of Kotliar, Lee 
and Read 26] ) . The distinction between these two situa- 
tions, together with a naive phenomenological extension 
of the large Ud problem, will be briefly discussed in Sec- 
tion 3.3. 

The model B focuses on the bonding band only, i.e. on 
an effective low-frequency description of the conduction 
electrons. So long as the underdoped regime is in ques- 
tion, the main low-frequency effects are expected to come 
from an opening of a correlation gap in the single-particle 
excitation spectrum. These effects are related here to the 
AF correlations and are described in the adiabatic way, 
by considering the influence of an effective single-particle 
interaction 

H COTI = ^[A(k) J Dj C(r £> k±Qff +h.c] (2) 
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on the real part of the electron self-energy, with Q = 
(n/a,ir/a) being the AF wave vector. Although the ef- 
fect of these correlations on the imaginary part of the 
electron self-energy is not treated explicitly, it can be 
represented by a phcnomcnological contribution to the 
damping energies which increases linearly with frequency 
(see a comment on the generalized Drude formula given 
below). The magnitude of A(k) is assumed to be real 
and small, as compared with the energies A p d and t p d- 

The structure of the Bloch functions and the Bloch 
energies of the model A is well known j2|| [2£j . On the 
other hand, the diagonalization of the model B (i.e. of 
the Hamiltonian + H co „) is straightforward, leading 
to the dimerized bonding band with the dispersions of 
the upper (L — A) and lower (L = S) subbands given by 

E A>s (k) = ±[E D (k) + EoQt)} 



±^ l [E D OL)-E R QL)]*+A*(tf, (3) 

where Ed_(\sl) = Eo(k ± Q). The effect of the perturba- 
tion Aeon- on the Bloch functions is given in the usual 
way, in terms of the auxiliary phase ip^ defined in Ap- 
pendix B. 



B. Coupling Hamiltonian 

According to Ref . jlE] , the coupling of the conduction 
electrons to the external electromagnetic fields, relevant 
to the in-plane optical conductivity analysis, is given by 
the coupling Hamiltonian 

H cxt = ~ £ K(q ± )J Q "'(k)4 +q ^i^ + h.c] 



LL'ka 
a 2 



2mc 2 



£K(q ± )(-) 7 ^(k;2)Ct +qi(T C kCT 



ku 



+h.c.]. 



(4) 



Here C is the index of the conduction band (C = D in the 
model A and C = A in the model B), and a S {x, y} is 
the photon polarization index. The J^ L (k) are the cou- 
pling functions in the first-order term (the current ver- 
tices) and (k; 2) is the coupling function in the bare 
second-order term (the bare Raman vertex). Further- 
more, A a (qj_) and A^(q±) are the Fourier transforms 
in space of the vector potential A a (r) and of its square 
A^(r), respectively, and c[± • a Q = 0. The coupling func- 
tions of the models A and B are shown explicitly in Ap- 
pendices A and B. 



C. Optical conductivity 

The doping dependence of the optical conductivity 
spectra of the I^-zSr^CuC^ single crystals has been 



systematically examined by Uchida et al. The un- 
derdoped crystals exhibit the behaviour which cannot be 
described as the response of a Drude metal. In particu- 
lar, two distinct threshold energies are found in the un- 
derdoped regime. The first one (MIR threshold energy) 
is strongly affected by the doping level (with the maxi- 
mum in the spectra at the energy liWp eak placed between 
0.1 and 0.3 eV). In the clean (underdamped) limit of 
the model B (Ao 3> S3, where Ao labels the magnitude 
of the parameter A(k) and £3 is the MIR damping en- 
ergy), this energy should be ascribed to the AF pseudo- 
gap (EaO^f) — Es(kp), with kp on the Fermi surface). 
In contrast to that, the aforementioned one-band models 
had success in explaining this structure in terms of the 
parameters U/2 [lj E2, J or E The second 

(charge-transfer) gap appears at 1.5—2 eV. In the model 
A, this energy scale is connected with the pd dimerization 
gap (-Ex(kp) — EoCkp), L = N, P), and, importantly, is 
not simply related to the MIR threshold energy. 



In the gauge invariant formalism |25l |42j , the optical 
conductivity of the underdoped high-T c cuprates can be 
thus shown in the form 



total 



i e 2 n? ff hu> 



fthj + i£i 



kja*(w, S 3 ) 



W Q ntc >, E 2 ) ~ ^KcoH - 1]. (5) 



The four quantities which enter into this expression rep- 
resent the effective number of conduction electrons per 
unit volume, n° s , the MIR polarizability, a"(ui, £3), the 
interband polarizability, a™ ter (a;, £2), and the contribu- 
tion of all other, on-site high-frequency optical processes, 

£a,oo(k-0- 

The £i are three phenomenological damping ener- 
gies. As pointed out above, the description of the low- 
frequency processes in Eq. JSJ) can be easily improved 
by extending the first term to include the frequency de- 
pendent corrections in the intraband damping energy 
(£1 — > £i(w) = H/r(u>), usually H/t(uj) = H/t + afiw) 
and the effective mass (which frequency dependence is 
forced by the causality principle, n" s /m — > n/m(u)). It 
should be recalled that the generalized Drude formula 
(a ^ 0, but c4 r (w) = and < ltcr (w) = 0) would hardly 
be extended to the underdoped regime because large 
damping energies are required to explain measured spec- 
tra Q. In the present analysis of the underdoped regime, 
where the MIR optical activity is associated with the ex- 
citations across the AF pseudogap, rather than with the 
strong intraband damping effects, such frequency correc- 
tions to £1 are presumably small, and will be disregarded 
in the quantitative analysis given in Section 3. 

In the model B, the first two terms in Eq. (0 are given 
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by 



ff = ^E^(k)[l-/i(k)], 



V 



(6) 



2^ 5 (k)[/ A (k)-l] 

(^ + i7j)2-^ s(k) - 



(7) 



The sum X)k* ' s restricted to the first Brillouin zone of 
the dimerized lattice, and 7^(k) is the static Raman 
vertex defined in Appendix B. Similarly, the model A 
gives rise to the interband polarizability of the form 



25 ^ E 2 LD (k) 
-2E LD (k)f D (k) 



-- y 

uj 2 V ^ 



(8) 



The Fermi-Dirac function [1 + e^ EL ^ ^] 1 is denoted 
by /i(k), and the energy difference £x(k) — Eli(\s) by 

Since the contribution of the on-site high-frequency 
processes is nearly independent of doping, and is small 
at energies below 3 eV, we take Im{e a!OC (u)} = and 
fie{^Q,oo( w )} ~ £oo m this energy region. (A better 
approximation for £ a ,oc{u) is given at the end of the 
article.) For further considerations it is appropriate 
to denote the first three contributions in Eq. © by 
a% Iudc (uj), o%(ij) and a™ teI (u), with the abbreviation 



_Drudc 



The corresponding macroscopic dielectric function 
reads as 



UJ 



(9) 



giving finally the three-component optical model with the 
minimal number of adjustable parameters: t p d, A p( j, A , 
Si and Eoo . (Actually, we take over the estimate of the 
ratio tpd/Apd from the EFG analysis |4l|, and reduce ad- 
ditionally the number of independent model parameters.) 



III. RESULTS AND DISCUSSION 

Let us remember some additional details about the 
optical anomalies of the underdoped and lightly-doped 
La2- a; Sr a ;Cu04 crystals, (i) A small transfer of the spec- 
tral weight across the charge-transfer gap is observed 
with decreasing doping in the 0.2 > S > 0.1 doping region 
(the derivative dflf DtcI /dS is negative but small; 0? nter is 
the interband spectral weight to be defined below), (ii) 
As mentioned above, this phenomenon is accompanied by 
the appearance of the MIR peak placed at hu>^ cak p=> 0.1 
eV. (iii) The additional decrease of doping, below <5 < 0.1, 
is characterized by a dramatic increase of the spectral 



weight above the charge-transfer gap (i.e. d£lf ntel /d6 < 
becomes large, resulting for 6 ~ in the complete disap- 
pearance of the intraband spectral weight), (iv) The MIR 
peak is also shifted, resulting in fiw p r cak ~ 0.15 eV and 
fiWpeak ~ 0.3 eV at 5 = 0.1 and 5 = 0.06, respectively. 

These features will be discussed now in the frame- 
work of the present three-component optical model. The 
considerations include three important questions. First, 
starting from the model parameters independent of dop- 
ing, we show how the total spectral weight is shared 
among the Drude, MIR and interband channels. Sec- 
ond, we illustrate the dependence of the low-frequency 
optical conductivity on the symmetry and magnitude of 
AF pseudogap. Finally, the validity of the model I© — © 
in the lightly doped regime is tested on comparing the ex- 
perimental spectra measured at too < 3 eV to the ones 
which correspond to the model parameters independent 
of doping. The characteristic energy scales used in the 
numerical calculations, A p d = 0.66 eV and t p d = 0.73 eV, 
are deduced from the estimate of the intracell hybridiza- 
tion t p d/ Apd ps 1.1 |4l| obtained on the basis of the intra- 
cell charge distributions measured in the La2-^Sr :c Cu04 
compounds by the EFG probes, and by assuming that 
the charge-transfer energy is roughly equal to 1.75 eV 
for the doping level S = 0.2. 



A. Conductivity sum rules 

There has been a lot of interest in the sum rules of the 
high-T c cuprates, in particular in the c-axis conductiv- 
ity sum rule and the in-plane Pines-Nozieres sum rule 
E3- Ei> both of which exhibit an indirect, relatively 
complicated dependence on the MIR threshold energy. 
Since it is not a trivial task in the three-band model to 
determine these sum rules, the related theoretical anal- 
yses have employed again the numerical treatments of 
the one-band models or suitable phenomenological mod- 
els [3 Eil ■ This contrasts with the direct dependence 
of the in-plane conductivity sum rule of the underdoped 
compounds on both the MIR and charge-transfer thresh- 
old energy, where it is possible to make a relatively simple 
estimation of the model parameters, in particular their 
dependence on the doping level. Let us now examine the 
latter issue in more detail. 

The experimental analyses present usually the total 
spectral weight of er^ ' 3,1 (uA in terms of the spectral func- 
tion N ctf (hoj) defined by gHU 



N eS (fru;) = -1 P 'da/ i?eK otal 
"o Jo 



(«')}• (io) 



The frequency fio = ^/iire 2 /(mVb) is a frequency-scale 
parameter, and Vq is the primitive cell volume. 

Alternatively, in the multi-component optical models, 
several auxiliary frequencies, together with the associ- 
ated effective numbers, can be defined by considering the 
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spectral weight contained by the i-th channel 



dw Re{a l a {io)} 



v nfn 2 . 



(11) 



(To simplify notation, we will refer to SI? as the spectral 
weight in the i-th channel. Obviously, this differs from 
the usual normalization of the spectral weight by a factor 
1/2.) The index i £ {total, intra, inter, Drude, ir}. No- 



tice that 

,eff 



eff 
t Drudci 



defined by this relationship, is identical 



to rtf in Eq. (JSJ. 

The conductivity sum rules of the present optical mod- 
els A and B are most appropriately described by 



"total, 



O 2 -4- O 2 

J 'intra, ' J 'inter, 0' 



O 2 — O 2 -J- o 2 

1 'intra — "Drudo^^'ir 



(12) 
(13) 



(the label in Eq. O refers to A(k) = 0). For the nar- 
row bands with the relaxation processes negligible (limit 
Si — * 0), the spectral function N c s(hij) has a step-like 
form, with the steps simply related to Eqs. I|12|) and Q13JI. 
However, in the typical experimental situations (where 
the bands are wide and the damping energies are not 
negligible) the relationship between these two represen- 
tations of the total spectral weight is more complicated 
(see Eqs. $Wl ). Nevertheless, a careful comparison be- 
tween the measured spectral functions N c ff(huj) and the 
analytic expressions for S7 2 can be used to extract the 
values of the parameters involved in the model JSJl , or in 
some more general optical model (see Ref. |44| . where 
the intraband spectral weight and the related effective 
number of electrons are estimated in various optimally 
doped cuprates and correlated with T c ). 

Shown in Fig. 1, the data A e ff(3eV) measured in 
La2_ 2; Sr :E Cu04 and the total spectral weight of the three- 
band model, 



O 2 

1 'total i 



.o = n o^£H7£? (k;2)Mk), (14) 



display the same doping dependence. At the same time, 
the doping dependence of the intraband spectral weight 
is given by 



O 2 

J 'intra, 



1 



(15) 



kcr 



The intra- and interband spectral weights are also shown 
in Fig.l. Interestingly, due to the gauge invariant form of 
cr^ otal (a;), none of these three spectral weights depends on 
the damping energies. (Here 7^f(k) is the static Raman 
vertex of the model A; see the effective mass theorem 

63)). 

It should also be recalled that this figure shows the 
results obtained in the simplest case, in which the pa- 
rameters A p d and t p d are assumed to be independent of 
doping. In this case, for the doping range relevant to 
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FIG. 1: The effective number of electrons in the model A 
as a function of the hole doping. The experimental points 
iVeff (1.5eV) (filled triangles), N e g (2eV) (open diamonds) and 
iV c ff (3eV) (filled diamonds) measured in La2-^Sr a ,Cu04 (Ref. 
4, with S = x assumed) are also shown. 



the hole-doped high-T c cuprates, < S < 0.3, the calcu- 
lated total spectral weight is shared between the intra- 
and interband channels nearly in equal proportions. 

Fit of the predicted spectral weights to the measured 
data is based on the relationships 



fi 2 

J 'total, 
O 2 

J 'intra, 



N cS (flU! 

N cS (huj- 



i ax)-\) ; 



(16) 



Here Sw max and huj m i n are the effective maximal and 
minimal energy in the interband processes, respectively. 
Although there is a large inaccuracy in the experimen- 
tal determination of fio> m i n and fia; ma x, it seems that 
0.2 < N c g(hxj ni i n ) < 0.25 and A^^^x) ss 0.5 almost 
in all La2CuC>4 based compounds 0, S (see, for ex- 
ample, Fig. 10 in Ref. 4]). In this respect, note that 
the MFA of the large Ud case (i.e. for the average bond 
energy involved both in the low- and high-frequency pro- 
cesses), together with the conclusions of the EFG anal- 
ysis, suggests the cutoff energy huj m i n 2 eV (Fig. 1), 
contrary to the experimental optical studies, which are 
consistent in the conclusion that fi^ m i n < 1.5 eV. In Sec- 
tion 3.3 we will turn to a more detailed discussion of this 
discrepancy. 

In the model B, one obtains that the spectral weight 
of the MIR processes is given by 



n 2 = 



\ ^ 



m2|J^ s (k)| 2 



k*cr 



EasM 



[1-/ A (k)], (17) 



and that of the Drude term by 



n 2 

5 'Drude 



V n 



off o 2 
"o- 



(18) 
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FIG. 2: The dependence of the Drude spectral weight on the 
symmetry of the AF pseudogap in the hole doped regime, for 
Ao = 50 meV. The filled diamonds represent the square of 
the measured plasma frequency Q shown in units of fio/^oo, 
as discussed in the text. The parameters are hQo w 3.8 eV 
and Eoc ~ 5. 



The resulting doping dependence of the Drude spectral 
weight is illustrated in Fig. 2, for three symmetries of 
the effective potential 



A(k) 
A(k) 

A(k) 



= A , 



-Ap(cosk • ai — cosk • a.2), 



-X^/i + i(cosk • ai - cosk • a 2 ) 2 , (19) 



corresponding, respectively, to the AF pseudogap of the 
s, d x 2_ y 2 and anisotropic-s symmetry. Again, t p d = 0.73 
eV is the average bond energy, and both spectral weights 
are independent of S^. 

The most noticeable feature of the Drude spectral 
weight is the occurrence of two different metallic be- 
haviors: (i) the electron-like behavior characterized by 
<?f2Drude/cW < 0, and (ii) the hole-like behavior where 
dflr>rude/d5 > 0. For A <C t p d, A p d, the critical doping 
<5o, at which the behavior of the conduction electrons will 
be changed, is only weakly dependent on the AF pseudo- 
gap symmetry. The critical doping associated with the 
curves in Fig. 2 agrees well with the results of the Hall 
effect measurements (where the sign r eversal of the Hall 
coefficient is found at 5o ~ 0.25 0,03). 

Finally, one can also notice that, when the doping level 
is reduced well below So, the effective number of charge 
carriers obtained by the expression © can be adequately 
described by the free-hole approximation. The effective 
number of holes is small with respect to the number of 
electrons in the conduction band, (1 — S)/Vq. Eq. (|18fl . 



together with the definition of the static Raman vertex, is 
the relevant expression to reconcile the concept of large 
Fermi surfaces, supported by the ARPES experiments, 
with the small effective number of holes, estimated from 
the Drude spectral weight or from the associated elec- 
tron plasma frequency (the measured plasma frequencies 
are indicated in Fig. 2, using a simplified single-band 



expression ^f) rudo ~ e c 



.(f^ + l/r 2 ) 



B. Low- frequency optical conductivity 

According to Fig. 2, the above description of the MIR 
correlations becomes ineffective for the magnitude Ao 
small in comparison with the energy difference between 
the Fermi energy and the van Hove energy, \i — e v h- 
In such circumstances, the two-component low-frequency 
optical conductivity, characterized by the distinct MIR 
threshold energy, will transform into a nearly single- 
component form with the increasing role of the frequency 
dependent corrections to the intraband damping energy 
[H 113 ■ Not surprisingly, these qualitative changes in 
the spectra are correlated with the sign reversal of the 
Hall coefficient |l8j . 

In order to express the dependence of the low- 
frequency optical conductivity on the model parameters 
more clearly, we show in Fig. 3 the normalized spectra 
for a few characteristic cases. First, it should be recalled 
that the intraband spectral weight is directly depen- 
dent on the electron group velocity J^ D (k)/e, namely 

^fntra.O = ^><ntra^0 and 



eff 
intra 



m 1 



iyL[C(kP[£D(k)-4 (20) 



As a result, the transfer of the spectral weight across the 
MIR gap will be strongly enhanced by the s component 
of the parameter A(k). However, a huge MIR peak has 
never been observed in I^-zSr^CuC^ or in I^CuC^+j; 
measured spectra (for example, compare the relative in- 
tensity of the MIR spectrum measured in La2Cu04.i2, 
the long-dashed line in Fig. 3(b), with the model predic- 
tion for A(k) = Ao, the long-dashed line in Fig. 3(a)). 
Therefore, it seems that A(k) is dominantly of the d x 2_ y 2 
character, with a small s component directly related to 
the velocity ratio 1/(1+A) measured in the ARPES in the 
nodal (k x = k y ) region of the Fermi surface pjil I3M l36l| . 
Second, in the present model, where the second-neighbor 
bond-energy t pp is set to zero, the anomalous shift of the 
MIR peak structure with decreasing doping can only be 
understood as a result of the doping dependent magni- 
tude Ao . The lower doping, the stronger effective pertur- 
bation has to be, in full agreement with the conclusion of 
a similar single-band analysis focused on the character- 
istic AF energy, given in Re f . [45f . and with the related 
specific heat measurements 46]. A significant deviation 
from such doping dependence of Ao is expected for t pp 
comparable to t p d, or for t pp Aq > ji — s v h- The 
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FIG. 3: The dependence of the normalized low- frequency op- 
tical conductivity on the AF pseudogap symmetry (a) and 
magnitude (b), and on the doping level (c), Ei = E3 = E. 
<7° c = ft^Drudc/ (47rE) is the DC conductivity. The exper- 
imental data measured in La2Cu04.i2 (long-dashed line in 
Fig. (b)) are given for comparison. In this material, 
ft<4 r eak « 90 meV and E « 30 meV at T w 200 K. 



C. Interband optical conductivity 

It is interesting to investigate finally how the expres- 
sion relates to the interband optical conductivity 
data measured in the lightly-doped and underdoped com- 
pounds. Note that both the U4 = case and the MFA 
of the large Ud case are unable to explain the anoma- 
lous doping dependence of the interband spectral weight 
(compare the model prediction for fif nter with N e g(3eV) 
-N cS (1.5eV) in Fig. 1). Also, when fitted with A pd = 
0.66 eV and t p d = 0.73 eV, the model © underestimates 
the interband spectral weight of the considered com- 
pounds nearly by a factor N c g (2eV) /N e g(1.5eV). The 
solution of the problem comes on noting that the MFA 
of the large Ud case treats the effects of strong local cor- 
relations on the bond energies in an average fashion. It 
is possible, however, to make a simple phenomenolog- 
ical extension of the present multiband optical model, 
and to improve the overall agreement between the model 
predictions and the measured data. This is due to the 
gauge invariant form of Eqs. ©— ©, which provides 
a successful separation of the twofold role of the inter- 
mediate states in the optical processes. For example, 
the copper-oxygen hopping processes participating in the 
low-frequency optical activity are accompanied by u> = 
and S2 = 0, resulting in the effective mass theorem (|A3I) . 
in full agreement with the longitudinal multiband re- 
sponse theory [2£j. Similarly, the high-frequency optical 
processes can be described in terms of the elastic Raman 
vertices (compare the interband polarizability with the 
difference 7 £f (k, u))-j££ (k) in Eq. JHJ|), with the DC 
conductivity behaving correctly in the metal-to-insulator 
phase transitions |2q . 

The parameters describing the low-frequency (high- 
frequency) optical processes in the strongly correlated 
systems are (are not) appreciably affected by the 
strong local correlations. The quantity that describes 
these effects is the probability of the photon absorp- 
tion/emission, in the low- frequency processes given by 
the current vertex J® D (k), in the high-frequency pro- 
cesses by J^ L {k). In the strongly correlated regime, 
these two vertices should be written in the form 



J, 



DDi 



(k) « &/Wa 00 and J^(k) « ^/t pd J Q ^(k), 
in the first approximation. Here t^ d is the bare bond en- 
ergy tpd i s the bond energy renormalized by the large 
U d , t p d is the average energy used in previous sections, 
and j£ L (k) are given by the expressions (| A 1 1) . Conse- 
quently, the optical conductivity of the large Ud case is 
more appropriately described by 



(21) 



difference in the shape of the Fermi surface between the 
model prediction for t pp = and the experimental find- 
ings can also be removed in this more realistic case. The 
low-frequency optical analysis of this interesting and ex- 
perimentally relevant situation will be given later. 



where 



-intra 



-inter 



M « (tZJtpdfa^u), 



(22) 
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energy (eV) 

FIG. 4: The real part of the optical conductivity for 5 = 0.06, 
the dj.2_ y 2-symmetry AF pseudogap with Ao = 0.125 eV, 
Ei = 50 meV, E 2 = 0.5 eV, E 3 = 0.25 eV, Ci = 0.2, & = 1.9 
and £3 = 0.8. The experimental data (open triangles) are 
from Ref. Q. 

or generally 

5iM « C^LM- (23) 

£1 , C2 7 and C3 a re three additional parameters attributed 
to the Drude, interband and MIR channels, respectively. 

To conclude, due to insufficient information regard- 
ing the low-frequency part of the experimental spectra 
in the lightly-doped crystals, the discrepancy between 
the model predictions with the average bond-energy t p d 
and the measured data is more clearly seen for the in- 
terband spectral weight. Its anomalous doping depen- 
dence can be understood as a result of a modification of 
the bond-energy renormalizations described by the factor 
C2 ~ (tpd/tpd) 2 in Eqs. J2J - (EH • A reasonable agree- 
ment of the data measured in the lightly-doped crystals, 
shown in Figs. 4 and 5, with the model Ij23(l is achieved 
for C2 « 2 (S = 0.06) and C2 « 3 (S = 0). (For com- 
parison, the value C2 ~ 1.3 is obtained by fitting the in- 
terband spectrum of Lai.gSrg^CuC^ to the model 
Not surprisingly, similar doping dependence of the bond- 
energy renormalizations is found in the previous slave- 
boson analyses [13,E3- For example, for the intracell 
hybridization &/A p d, — 1.1, one obtains (tp d /tp d ) 2 w 3.5 
|4?|. in a qualitative agreement with the present conclu- 
sions. 

Obviously, for the low-frequency part of the spectra, 
this is an oversimplified view, because the effects of the 
strong local correlations on these processes are unlikely 
to be described simply by multiplying the MFA spectra 
by the factors £i and £3. To clarify this essential ques- 
tion, further experiments in different high-T c families are 
desired as a function of the doping level and with a better 



12 3 

energy (eV) 

FIG. 5: Main figure: The real part of the optical conduc- 
tivity at various £2 for 8 ~ and fi = £3 = 0, compared 
with the experimental data (open triangles) |4|. Inset of fig- 
ure: The on-site contributions to the optical conductivity 
aj/(47r)Jm{e Qi0O (aj)} (i.e. the three-band-model contributions 
subtracted from the experimental spectra) for S — 0.06 (solid 
line) and 5 = (dotted line). 

resolution of the low-frequency part of the spectra. The- 
oretically, one needs more accurate investigation of the 
relations between the low-frequency optical conductivity 
and the Raman and ARPES spectra. 

Finally, notice that the difference between the exper- 
imental spectra and the conductivity cr^ otal (a->) at en- 
ergies higher then 2 eV can be ascribed to the on- 
site interband optical conductivity i?e{<7° n-slte (w)} = 
w/(47r)7m{£ Qi00 (w)}. The on-site spectra deduced from 
the experimental data and the three-band model spectra 
(solid lines in Fig. 4 and 5) are shown in the inset of 
Fig. 5, for the cutoff energy huj ma ,x ~ 3 eV, and for 5 = 
and 0.06. The nearly identical frequency dependence of 
Re{a° n ~ slte (u>)} for these two doping levels gives an ad- 
ditional, although rough support to the picture of the 
interband optical conductivity discussed above. 



IV. CONCLUSION 

In this article, the influence of two dimerization po- 
tentials on the optical response of the valence electrons 
in the high-T c cuprates has been studied, starting from 
the model parameters independent of doping. It is found 
that for the nearly half filled conduction band the pd 
dimerization potential leads to the spectral weights of 
the intra- and interband channels nearly equal to each 
other. It is also shown that at optimum doping the 
effective correlation potential related to the AF pseu- 
dogap processes results in the non-Drude nearly single- 
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component low-frequency optical conductivity. It is ar- 
gued that, due to the decrease in the effective number 
of conduction electrons, the low-frequency conductivity 
transforms naturally from a single-component into a two- 
component structure with decreasing doping. Both of 
these results are found to be in accordance with experi- 
mental observation in the lightly-doped and underdoped 
La2Cu04 based compounds. In the present model (where 
the bond energy t pp is set to zero), the pronounced dop- 
ing dependence of the MIR peak structure is attributed 
to the doping dependent magnitude of the characteris- 
tic AF energy. Finally, the substantial increase of the 
interband spectral weight at small doping levels is inter- 
preted in the framework of the conductivity sum rules 
as an indirect evidence on the strong renormalizations of 
the parameters describing the low-frequency physics. 




FIG. 6: The usual diagrammatic representation of the elastic 
Raman vertex, L ^ C. In the static limit of the model A (the 
index C = D), the summation gives the expression lA3t . 
i.e. the inverse k-dependent effective mass. The wavy and 
solid lines represent, respectively, the photon and electron 
Green functions. Notice the minus sign in front of the Raman 
vertices which reflects the fact that electron-like conduction 
band is shown in the hole picture. 
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APPENDIX A: COUPLING FUNCTIONS IN THE 
MODEL A 

By using the auxiliary functions Uk, «k and ik defined 
in Ref. [29j, one obtains that the coupling functions rel- 
evant to the case of the nearly half filled bonding band 
are 



Jf p (k) 



eat pd 2M k i> k 



eat 



pd M k 



tk 



sin k ■ 



• sin k ■ 



eat 



J 



and 



The 



DN (k) 

?"00 - - 

nass scale 



P d 2wk . 1 , 1 , 

— sin -k • a 2 cos -k • ai, 

H t k 2 2 

eat 2 pd 2u k . 1 1 
— p— sm -k • ai cos -k • a 2 , (Al) 

tl tk z z 



m 



1, 



sin -k • 

fTlxx 2 

has the usual form 



(A2) 



H 2 A pd /(2a 2 t 2 pd ), with |ai| = |a 2 | = a. 

By comparing two expressions for the dimensionless 
static Raman tensor obtained in the static limit of the 
longitudinal [2£j and transverse response of the model A, 
we can verify the effective mass theorem 

to d 2 E D {k) 



= (-) 



h 2 dk 2 



( k ; 2 ) + 3 E 

L=P.N 



2|Jr(k)h 
E LD (k) 



(A3) 



The related inverse k-dependent effective mass reads as 

(l/TO) 7 £?(k). 

Beyond the static approximation, the effective second- 
order term in the electron-photon coupling comes from 
the contributions shown in Fig. 6. The elastic Raman 
vertex is the vertex function in such effective processes. 
The gauge invariant form of the elastic Raman vertex 
reads as 



t£?00 



m 

72 



E 

L=P,N 



E 2 LD (k) 



2E LD (k) 



(Hu) + irf) 2 



E 2 LD (k)- 



(A4) 



It should be noted that the usual form of the elastic Ra- 
man vertex [3(], EE OH comes on disregarding the factor 
(hw) 2 /El D (k) in the expression {53} . 



APPENDIX B: 



COUPLING FUNCTIONS IN THE 
MODEL B 



The influence of the perturbation H COTT on the Bloch 
functions can be shown in terms of the auxiliary phase 
which satisfies the relations 



COS (y5k = 

sin <y9k = 



E D (k)-E R (k) 

EasM 
2A(k) 



EasW 



(Bl) 



The resulting static Raman vertex, the elastic Raman 
vertex, and the current vertices relevant to the < <5 < 1 
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doping range are given by and 

7^(k) = 7^ (k) cos 2 <f + 7 £^(k) sin 2 f J^W = 00 -S 2 ^ + J^(k) sin 2 If 



2|J^(k)| 2 « J^(k)cos^ k , 



e 2 ^s(k) ^ ^ J^( k ) = [J^( k) - j££ (k) ] sin f cos ^ 

7 r(k)co S ^ k -g 2| ^ g S |2 , (B2) « Jr(k)sin^ k . " ~(B4) 



e 2 S AS (k) 
) 2 |J a AS (k)| 2 

£2 s (k) £) band^of the model A, and 7^ (k; 2) « 7fcP( k ) cos ip k 



AA AA/-I \ m (Hlo) 2 \J^ s (h)\ 2 Again, the labels D and k refer to the k±Q states in the 

l aa ( k ,^J = laa ( k ) + -5 r2 n,\ D hand of the model A. and V^fk: 9.) M 



2SAs(k) 



111 



(B3) 



Eq. gj. 
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